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Abstract 

The Titius-Bode law for planetary distances is reviewed. A model describing the basic 
features of this rule in the "quantum-like" language of a wave equation is proposed. Some 
considerations about the 't Hooft idea on the quantum behaviour of deterministic systems with 
dissipation are discussed. 
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1 Introduction and motivation 

In the last years many papers, in particular those of 't Hooft [T], but also of several other authors 
[2J [3J U], have described a possible way out from some of the unpleasant (or, to cite Feynman, 
"peculiar") features of Quantum Mechanics. 

The idea proposed in [T] is, in a nutshell, that a classical deterministic theory supplemented with 
a dissipation mechanism (information loss) should produce at larger scales a quantum mechanical 
behaviour. The emergent model, after that the dissipation mechanism has been implemented, can 
be described with the usual tools of Quantum Mechanics, namely, states evolving in Hilbert space, 
unitary matrix, wave equation, etc. 

Originally motivated (also) by the challenge of Quantum Gravity, the initial accent on the Planck 
scale collocation of the deterministic underlying theory has been more recently weakened in favour 
of the description of several classes of deterministic classical systems, where quantum behaviour 
emerges from classical dynamics constrained by dissipation mechanism(s) (see [2H31EI)- So far, no 
explicit scale has been implemented in such models. 

On the other hand, it is quite clear that dissipative/chaotic processes (of energy, angular momen- 
tum, etc.) have acted, and act, also on scales very far from the Planck scale. For example, in the 
proto-planetary nebulae there are many dissipative and chaotic processes due to friction, viscosity, 
gravitational resonances, etc. Therefore, it is worth asking if these dissipative/chaotic processes 
could produce, on different scales, results analog to those predicted by the 't Hooft theory, namely 
some sort of quantum-like behaviour. Surprisingly enough, Nature seems to play just in this way. 
The Titius-Bode law of planetary distances seems to bring evidence along this direction. In analogy 
with the 't Hooft idea, the chaotic dissipation occurred during the formation period of a planetary 
system could have produced "the apparent quantization of the orbits which resembles the structures 
seen in the real world" [5]. In the present paper, an "effective" quantum-like kinematics (in the 
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form of a deformed Schrodinger wave equation) is developed and applied to the description of plan- 
etary systems. In doing this, we obtain also a quite precise prediction of the internal radius of the 
rings surrounding Saturn, Jupiter and Uranus. In the end, some qualitative considerations about 
permitted orbits and dissipative chaotic mechanism(s) are discussed. 

2 Titius-Bode law 

The Titius-Bode law is the most famous among the relationships about the spacing of planetary 
systems. It is an empirical rule which gives, under the hypothesis of circular orbits, the distances 
of the planets from the Sun as a function of a single parameter, i.e. an integer n. There are several 
versions of the law. The eldest form is perhaps the following 

r(n) = 0.4 + 0.3 ■ 2" (2.1) 

where r(n) is given in Astronomical Unit (1 A.U. ~ 150 • 10 6 km). For n = — oo, 0, 1, 2, ... this rule 
gives the distances, respectively, of Mercury, Venus, Earth, Mars, etc., including the asteroids belt 
(actually, Cerere was discovered following the indications of this law) and Uranus, which at the 
moment of the first formulation of the law (1766-1772) had not yet been discovered (see the book 
of Nieto [6] for history, explicit formulations, theories and extensive comments). 
In this original formulation, the law was not able to account for the distance of Neptune and Pluto. 
The more recent versions of the rule, elaborated during the XX century by Blagg (1913) and 
Richardson (1943), are able to describe not only the planetary distances within the solar system, 
including planets like Neptune and Pluto, but also can be successfully applied to the systems of 
satellites orbiting Jupiter, Saturn and Uranus. The agreement between the predicted and the 
observed distances of the various satellites from the central body is really good, of the order of 
a few percents, as can be checked in the tables of Nieto. Of course, criticisms have been risen 
even against an effective physical meaning of the law (see, for example, [7]). However, in [HJ, it is 
also convincingly argued that the agreement with the observations cannot be safely considered as 
a mere statistical chance. This modern version of the Titius-Bode-Richardson rule seems then to 
have a quite firm status, certainly not yet comparable to that of the Balmer formula for the atomic 
spectra, but at least similar. 

The main feature of these versions of the Titius-Bode law is that the rule can be expressed, if we 
neglect second order corrections, by an exponential relation as 

r = ae 2Xn , (2.2) 

where the factor 2 is introduced for convenience reasons and n = 1, 2, 3, 

For the Solar System we have 

2A = 0.53707, e 2A ~ 1.7110, 

a = 0.21363 A.U. 

The amazing thing found by Blagg was that the geometric progression ratio e 2A is roughly the same 
for the Solar System and for the satellite systems of Jupiter (e 2A ~ 1.7277), Saturn (e 2A ~ 1.5967), 
and Uranus (e 2A ~ 1.4662). The parameter A is dimensionless, its value is taken from the observed 
data, for any of the systems considered (planetary or satellite systems) . The parameter a has the 
dimension of a length and it is linked to the radius of the first orbit of the system considered (in 
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fact r(l) = ae 2X ). 

Many theories have been developed during the last 240 years in order to explain, or describe, the 
Titius-Bode law. 

There have been dynamical models connected with the theory of the origin of the solar system, 
electromagnetic theories, gravitational theories, nebular theories. They all can be found in literature 
(see, for example, |9J) and they have been excellently reviewed in the book [6]. We remind that 
also the idea of using a Schrodinger-type equation in order to give account of the law (|2.2j) is not 
new: during the years many authors have suggested various approaches in this direction (see, for 
example, an incomplete list of papers in [10J jj. Also the present paper aims to develop a quantum- 
like model of the T-B rule, but this will be done through a deformation of the usual Schrodinger 
equation. 



3 Bohr model of the hydrogen atom 

In this section we remind briefly the Bohr model for the hydrogen atom, just to set the notation. 
The electron orbits are supposed circular (this will be held for planetary orbits also). The two main 
equations are the equation for the force (i.e. the equation of motion) 

„2 

m- = % (3.1) 

where "m" and "e" are the mass and charge of the electron; and the quantization condition on the 
(^-component) of the angular momentum 

mvr = nh n = 1, 2, 3, (3-2) 

In the Bohr model, all the orbits belong to the same plane, and this is also taken for true in the 
planetary models. From the two equations above, one easily derives 

h 2 „ 

n 



me 2 



(3.3) 

hn 
mr 

The first equation is the law of electron distance from the nucleus in the Bohr model. With this 
law, from the classical expression for the total energy we get the energy spectrum of the bound 
orbits 

1 9 e 2 e 2 me 4 1 

E = -mv 2 = = -r- -jr (3.4) 

2 r 2r 2h 2 n 2 v ; 

Now we shall try to apply analog ideas to the solar and planetary systems. 



1 However, in almost any of the papers of Ref. [TO], the Newtonian potential 1/r is usually assumed together with 
a Bohr-like quantization condition, and these two things together yield (as it is well known) a law for the orbital radii 
as r ~ n 2 (not as r ~ e 2A "). This forces the authors either to use non consecutive integers to accomodate the planet 
distances, or to consider the solar system divided in two subfamilies (i.e. terrestrial planets and gigantic planets). 
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4 Model a la Bohr for a planetary system 



In this section we introduce a model for the "quantization" of a planetary system. The model 
acquires its discrete, or "quantum", properties from a modification of the Bohr quantization rule 
for the angular momentum. The equations here proposed, for a generic planet of mass m, orbiting 
a central body of mass M, are 

v 2 GMm 
m — = „ — 

(4.1) 

— = vr = s e Xn 
m 

where n = 1,2, 3, s is a constant, and J is the angular momentum of the planet. Some comments 
are immediately required: 

• Because of the principle of equivalence, the masses "m" on the LHS and on the RHS of the 
first eq. (|4.1|) cancel out each other. 



• The constant s in the RHS of the second eq. (|4.1j) has the dimensions of an action per unit 
mass. It plays the role of h and it must be understood as an action typical of the planetary 
system under consideration. It is not possible to use h itself, because this would fix the wrong 
initial radius in the Titius-Bode law, that is the constant a in r = ae 2Xn . 

• The constant A in (j4. If) is the phenomenological parameter obtained from the observation 
(2A = 0.53707 for the Sun, 2A = 0.54677 for Jupiter, 2A = 0.46794 for Saturn, 2A = 0.38271 
for Uranus). 

• In the second of the eqs. (14. lj) . we quantize the angular momentum per unit mass (i.e. the 
velocity field and not the momentum field). This a consequence of the principle of equivalence. 
If we did not do so, we would obtain a law for r(n) where the scale of distance changes from 
a planet to another, as the planetary masses change from a planet to another. 

From (14.11) one immediately gets 

o p An „2 
se s „2An 



which is the Titius-Bode law if we identify a = s 2 /(GM). 

We can also compute the energy spectrum for the i-th planet from (14.2 



1 2 GMuu GMxiH (GM\ 2 m* 

where n = 1,2,3, ... . As we see, the energy of the i-th planet is not properly quantized by itself. 
This is because the mass im changes in general with the planet, and this implies different sets of 
energy levels for different planets. Instead, the energy per unit mass 

g . = ^ = _(GMyj_ = _l_GMJ_ 

m; \ s J 2e 2Xn 2 a e 2Xn y ' ; 
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is exactly quantized, i.e. it is a quantity which depends on the orbit n only (apart from the general 
constants G, M, s), and not on the mass of the planet belonging to that orbit. Therefore, the 
energy-per-unit-mass levels are valid for the whole set of orbits of a specific planetary system. 

Also here some comments are needed, in order to complete the explanation given before. 

• The constant s can be computed in terms of the mass M of the central body and of the 
parameter a (remind that the radius of the first orbit is r\ = ae 2X ) 

s = \jGMa (4.5) 



• The constant s is not the same for all the planetary systems (Sun, Jupiter, Saturn, Uranus). 
In fact, if it were so, this would imply that the parameter a = s 2 /(GM) should be in inverse 
proportion to the mass M of the central body, which is not true. Therefore the constant s is 
not universal, like h, but it depends on the planetary system under consideration (it changes 
from planetary system to planetary system). 

• In the Titius-Bode law the mass of the single planet (or satellite) plays no role, essentially 
because of the principle of equivalence. A naive quantization a la Bohr, with a quantum 
condition mvr = se Xn , would yield a system of orbits with radii 

(and corresponding energy levels). In such a system, every specific i-th planet, with mass 
nij, would have its own system of permitted orbits, and the coefficient a = s 2 /(GMmj) of 
the T-B law would depend on the mass of the specific planet under consideration. All this is 
contrary to the observed evidence. In a planetary system the coefficient a of the T-B law is 
the same for all the planets belonging to such system. (Of course, a changes from planetary 
system to planetary system). Agreement with the observed evidence can be obtained if we 
"quantize" the angular momentum per unit mass, J/m = vr = se Xn . 

The same considerations hold for the energy levels. The set of energy levels of the i-th planet is 
given by 

*w = - ass?- («) 

This implies that every planet has its own set of energy levels. Therefore, in order to have a 
quantity linked to the orbit only, and not to the mass of the planet belonging to that orbit, we 
should consider the energy per unit mass, £ = Eijm^. 

The difference between a planetary system and a muonic/electronic atom is in this respect striking. 
In a muonic atom the levels belonging to the muon can be occupied by the muon only; the electron 
has different levels. In a planetary system the orbit is singled out by the orbital velocity only. 
Jupiter can be put in the orbit of the Earth, if to Jupiter is given the same speed of the Earth. 
This is a simple consequence of the principle of equivalence. 
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5 Wave equation for a planetary system 



The Bohr-like model introduced in section [J] allows us to look for the corresponding Schrodinger-like 
equation. On postulating the absence of interference phenomena among the probability amplitudes 
of the planets along their orbits, and of course the absence of quantum jumps (i.e. transitions) 
among the orbits, we shall now construct a wave equation, and we shall see what it predicts for 
the energy spectrum and for the exponential positions of the planets, when the usual probabilistic 
interpretation of the wave function ip(x,i) is adopted. 

We suppose, hence, to associate to a planetary system a scalar field ip(x,t), the so called wave 
function. We shall see that the wave function does not give us information on the behaviour of 
the single planet, but rather on the global structure (orbits, energy levels, rings, etc.) of the whole 
planetary system (about this, see the papers of Reinisch and Nottale in [10]). 



On comparing the Bohr quantization condition mvr = nh with the condition given in (|4.ip . 
vr = se Xn , we see that the most straightforward correspondence is 



h 

s < ► — 

m 



or 



h < > sm. 



(5.1) 



This correspondence allows us to write down immediately a wave equation for the stationary states. 
In fact, from the Schrodinger eigenvalue equation we can write 



and defining 



£ 



V(r) 



we arrive at 



V 2 



E 



U(r) 



tp = Eip 



(5.2) 



— = energy per unit mass 
m 

u 

— = potential energy per unit mass 
m 



(5.3) 



-V 2 + V{r) 



i> = Sip, 



(5.4) 



which will be the starting point for the construction of our wave equation for planetary systems. 
We note again that the quantity correctly quantized is £ , the energy per unit mass, and not the 
energy itself. As already said for the problem a la Bohr, this is a direct consequence of the principle 
of equivalence. 



6 Wave equation for the Titius-Bode problem 

In this section we look for a wave equation corresponding to the model (14. ip . We also take into 
account the observed fact that all the planetary orbits lie (more or less) in the same plane by 
considering a 2-dimensional Schrodinger-like equation, i.e. we write the previous wave equation 
in a plane (on this, see also M. de Oliveira Neto et al. in Ref. pH]). This choice will suggest to 
us how to modify the standard Schrodinger equation in order to accommodate for the particular 
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quantization condition (|4.ip on the angular momentum, namely vr = se Xn . 
The equation (|5.4p can be written in operatorial form 



H M tp = £tp 

where H M is the Hamiltonian per unit mass 

Hu = \ + V(r) 

with the association p <-> — is V. 

The Hamiltonian H M , written in two dimensions and in planar polar coordinates, reads 



(6.1) 



(6.2) 



(6.3) 



With the usual associations of ordinary quantum mechanics, in planar polar coordinates, 

„2 ■> 1 9 ( d 



Vr 
P<p 



2 1 d 
r dr I dr 
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the Hamiltonian (|6.3p reads 



s 

2^2 



-18- 



r d_ fd_\ d 2 
dr \ dr J d(j) 2 



+ V(r) 



(6.4) 



(6.5) 



The operator which accounts for the whole angular momentum is p^ (in fact, the motion takes place 
in the plane (x,y), that is, the only non zero component of the angular momentum is J z ). On the 
Hilbert space £ 2 ([0, 2tt]), the operator^ is hermitian (self-adjoint) with eigenfunctions e %m ^/^/2~K 
and eigenvalues ms, m £ Z {j>^e %m ^ = ms e 1 '" 1 ^). This corresponds to the quantum condition (a la 
Bohr) vr = ms (see (|3.2p ). Instead, in order to accomodate as eigenvalues the numbers se Xm , and, 
of course, to do not spoil or loose the eigenfunctions e tm ^ and their Hilbert space £ 2 ([0, 2ir\), we 
define the operator 



i e mX e im< ^ 



(6.6) 



for every e tm ^ € £ 2 ([0, 2tt]), m € Z, where A is the phenomenological parameter given in 
(2A ~ 0.53707 for the Solar System). 

Note that such a definition is directly dictated by the model (|4.ip . if one wants to reproduce the 
correct eigenvalues of the angular momentum-per- unit-mass (as we learnt them from (|4.ip ). The 
operator V v can be defined to be linear (see Appendix), and being defined on an orthonormal 
basis of £ 2 ([0,27r]), is therefore well defined on all £ 2 ([0,2-7r]). In the Appendix we show that the 
operator V v verifies 



ie 



Thus, we change the usual association 
definition of p^ we have 



-isdsh into the new one 



-isV^e 



im<f> 



se 



se \m e im4 



(6.7) 

isVtp. With this new 



(6.8) 
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The operator p^ so defined can be proved to be self-adjoint (see Appendix). 

Evidently, being the operator Vtp constructed by A, it will acquire its full physical meaning in a 
framework where the constant A is derived from a chaotic/dissipative mechanism at work in the 
proto-planetary nebula. In this direction, contacts are expected, for example, with the formulations 
proposed by Laskar, or von Weizacker [9], or Willerding [TT] . 
The new Hamiltonian operator per unit mass reads 



H - S " 
2r z 



9 ( 9 > r, 



+ V(r). (6.9) 



dr \ dr / 

The wave function in these planar polar coordinates can be written as 

ip(r,(f)) = R(r)<5>(<p) (6.10) 

where the separation of variables is supposed, in order to proceed towards a solution of the 
Schrodinger-like equation. The normalization condition is 

/•OO /*27T 

\ 2 d 2 x = R 2 (r)rdr- |$(0)| 2 # = 1. (6.11) 
Jo Jo 

The time-independent Schrodinger-like equation is H M ip = Since it is a 2-dimensional equation 
and it involves a deformed angular momentum operator, we proceed now to work out the solutions 
explicitly. With the position ip(r, 4>) = R{r)Q{<ft) the equation becomes 

r 3 ( dR\ 1^2^ 2r2 ,n T ,/ ^ 



and separating the variables we get the radial and angular equations 

(6.12) 



1 d f dR\ n n , 2 
dr \ dr J 



Z—\ r —)-^2 R +-( £ - V{r))R = 



where we suppose \i S 



7 The radial equation 

The radial equation is quite similar to the standard one of the hydrogen atom theory, apart of 
course for the radial part of the Laplacian, which is here 2-dimensional. For its solution we shall 
therefore use standard techniques (see e.g. [13]). 

Let's now look for the asymptotic behaviour of R(r). We ask R(r) to be finite everywhere including 
r = 0. Under the hypothesis 

HmV(r)r 2 =0, (7.1) 

1 — *0 

which is fulfilled by V(r) = —GM/r, the radial equation for r — ► becomes 

4 -"»="■ (72 » 
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We seek R(r) in the form of a power series and we retain only the first term for small r. That is, 
we put R{r) = kr l for r — > 0. Substituting this in the equation (|7.2p we find 

t 2 = p. (7.3) 

We want R(r) real, therefore t must be real, and fj, > 0. So we have two roots 

h = -y/fj- 

t 2 = +VJ1- (7.4) 

But t\ < 0, hence r* 1 — ► oo for r — > 0. So t\ does not yields a R(r) finite near the origin, and must 
be discarded. The only acceptable solution is ti = +y / 7i > 0. Therefore we put R(r) ~ kr t2 for 
r — > 0. For the Newtonian potential V(r) = —GM/r the radial equation reads 

or z r or r z s z \ r y 

We choose as natural units for mass, length, and energy (per unit mass), respectively, 

s 2 C?_M 2 
GM' s : 
so that the radial equation can be rewritten as 

d 2 R ldR u n n / 1\ „ 

_ + __-iii2 + 2 £ + - LR = 0. 7.6 

or z r or r z \ r J 

To study the discrete spectrum (bound orbits, £ < 0), we introduce, in place of £ and r, the 
variables 

1 2 

n = and p = -r = 2\/-2£r (7.7) 

with £ < 0, n > 0, and their inverse relations 



In 

£ = -^ r= 2 p - (7 " 



Equation (|7.6p then becomes 



d 2 R 18R ( 1 n p\ „ . . 

We already know that ~ /?* 2 for /) — > 0. If now we take p — ► +oo, then eq. (tXSD reads! 

fl2n i 

whose solutions are e ±P//2 . We want R{p) — > for p — > +oo, therefore we must choose the second, 

i?(p) ~ e _p/2 for p -> +oo. 

If now we make the substitution 

R(p) = pt* e~ p/2 w{p) (7.11) 



2 Since the normalization condition (|6.11[l implies R(p) — > 0, then R'(p) — ► 0, for p — > +oo. 
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then the eq. (|7.9p becomes 

pw" + (2 t 2 + 1 - p) w' + (n - t 2 - -) w = 0. 

We look for a solution of this equation which diverges at infinity no more rapidly than a finite power 
of p, while for /)->0we should have w — ► u; Q finite. Such a solution is the confluent hypergeometric 
function (see Appendix) 

= F(Q, 7 ,p) = F(i 2 + X - - n, 2t 2 + 1, p). (7.12) 

In particular, it behaves as a polynomial (w —* pP for p — » +oo) only if a = —N with N > 
integer. Thus 

t 2 + ^-n = -iV n = t 2 + ^ + A r , (7.13) 

where N = 0,1,2,3,.... 

8 The angular equation and the spectrum of 8 

We need now to know what is t 2 = yfji. To this aim, we solve the angular equation 

T>l<S> = -n$ p>0. (8.1) 



We take $ G £ 2 ([0,2tt]) and we know that an orthonormal basis in £ 2 ([0, 2ir\) is e !ra */\/2i, m £ Z. 
From the definition (|6.6p we immediately get 



which means that 

• the eigenvectors of "P 2 are u m = e im< ^ / V2tt; 

• the eigenvalues of V 2 are {— e 2Am } mg ^, and therefore yu m = e 2Am , m = 0, ±1, ±2, .... 
Hence we have 

t 2 = ^/p = e Xm m = 0,±l,±2, ... (8.3) 

and then 

n = t 2 + ^ + N = e Am + I + iV (8.4) 
with N = 0,1,2,3,..., m = 0, ±1, ±2, .... 

Finally, we are able to write down the spectrum of the energy (per unit mass) which in ordinary 
units reads 

1 / GM \ 2 1 1 GM 



£ 2\ s ) (e Xm + l/2 + N) 2 2 a (e Xm + 1/2 + N) 2 ^'^ 

with N = 0,1, 2, 3, m = 0, ±1, ±2, ... . 
A couple of considerations are now in order: 
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• For ./V fixed, and large positive m, the levels behave like 

1 / GM\ 2 1 , , 

i.e. we recover the formula (|4.4p obtained from the model a la Bohr. 

• Clearly, the best fit between the relations (|8.5|) and (|4.4|) (the relation (|4.4|) is directly given, 
in practice, by the observation) is obtained for N = 0. The sequence N = and m = 1, 2, 3, ... 
is called the principal sequence. The formula (|8.5p seems to suggest the existence of other 
sequences. In particular, we shall see that the sequence 

iV = and m = 0, -1, -2, -3, ... . (8.7) 

can be interpreted a system of rings (see section ITU]) . 



9 Mean value of r 

In order to complete our analysis of the solution of eq. (|6.12p , we want to compute the mean values 
taken by the variable r in various eigenstates. We are particularly interested in the eigenstates 
of the principal sequence, N = 0, m = 1,2,3,.., which is the one matching the observational 
data in the closest way. For tp = i?(r) <£(</>), the normalization condition (|6.1ip holds, and the 
= e trni t > /\/2ir are already normalized to unity. Therefore we are left with the condition on 
R(r) only 



1 = / drrR(r) z . (9.1) 
J o 

First, we find the correct normalization constant for R{p) by translating the normalization condition 
for R(r) into the one for R(p). Writing 

R(p) = A p t2 e~ p/2 w(p) , (9.2) 

we have (reminding p = 2r/n) 

r°o 2 rco 

1= drrR(rf = --A 2 dp p 2t2+1 e~ p ' [w{p)f . (9.3) 
Jo 4 Jo 

Since we are mainly interested in the mean values of r for the eigenstates belonging to the principal 
sequence, (N = 0, m = 1, 2, 3, ..), we set N = n — ti — 1/2 = 0. This means n = ti + 1/2 and 

w(p) = F(0, 2*2 + 1, p) = 1 (9.4) 

This simplify very much the calculation of the integral. Finally we have 

A 2 = „ - 2 (9.5) 

n 3 r(2t 2 + l) 

where T is the Euler' T-function, and we used 2n = 2ti + 1 and Y{x + 1) = xT(x). 
We can now compute the mean value of r: 

„3 poo 

drr 2 R{rf = —A 2 / dpp 2t ' 2+2 e -P [w(p)] 2 . (9.6) 
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Again, we are interested in the principal sequence. Setting N = and w(p) = 1, and using for A 
the value just obtained in (|9.5p . we have 



o /*oo o 

r = —A 2 d pp^ +2 e-P = —A 2 T(2t 2 + 2,)=n 2 + -. (9.7) 
8 J 8 2 

Since N = 0, then n = i 2 + 1/2 = e Am + 1/2. Therefore 

r = n 2 + ^ = (e Am + 1/2) 2 + ^(e Am + 1/2) (9.8) 
Restoring the ordinary units, we have as the mean value of r, for large and positive m, 

r ~ -^-e 2Am = ae 2Xm (9.9) 
which agrees with the Bohr model developed in section [H 

We see here that the agreement with the observed T-B rule holds only for large m. The error 
is particularly relevant for the internal planets, a feature unluckly shared with all the "quantum" 
models of the T-B rule. In some of them (e.g., those of Reinisch, Nottale, Oliveira) non consecutive 
integers have been chosen in order to recover contact with the observation. Others models (Al- 
beverio) also cite a necessary breakdown of the "quantum" T-B rule near the sun, because effects 
different from the gravitational ones become important. We tend to support this last explanation. 



10 Description of the rings 

It is tempting to speculate on other possible sequences, in particular the one N = and m 
0, —1, —2, —3, ... . First, we note that for N = the mean value r is still given by the relation ([9 
or in ordinary units 



(10.1) 



(e Am + l/2) 2 + i(e Am + l/2) 

where a = s 2 /{GM). A we said, the formula (jlO.ip matches the phenomenological formula (|2.2|) 
only for large, positive m. And this could be expected as due to the "quantum" origin (i.e. from a 
wave equation) of the formula (jlO.ip . But the interesting feature of the equation (jlO.ip is that it 
can be now considered also for m = 0, —1, —2, —3, .... In the limit m — ► — oo we have for (jlO.ip 

f — ^a, (10.2) 

contrary to the limit given by the phenomenological rule (|2.2p . which predicts r — * for n — > — oo. 
The sequence m = 0,-1,-2,-3,... would in such a way correspond to a system of permitted 
concentric orbits, accumulating on the limit orbit r = a/2: clearly, a system of rings. We can 
check the predictive ability of the relation (|10.2p using the planets equipped with a system of rings: 
Jupiter, Saturn, Uranus. The procedure is the following. 

First, we calculate the parameter a from the phenomenological Titius-Bode rule, using the 
radius r\ of the first satellite orbit: 

a = %. (10.3) 
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Rin—ring — n Q — n_9\ ' (1^.4) 



Then we calculate the radius that the inner ring should have: 

1 r\ 

Finally, we compare the above relation with the observational data, using the reverse form 

2e Ri n — r i n g = T\. (10.5) 

With the data given, for example, in [12] we can write, for: 



Jupiter: 



2 e 2A Ri n -ring = 2 • 1.7277 • 122 500 = 423 280 km 



This value agrees almost perfectly (within an error of less than 1 %) with the radius of the orbit 
of Io, which is 421600 km. As first ring, we considered the ring "Main". Besides, we considered Io 
as the first satellite (as regard the distance from Jupiter), instead of Metis, Adrastea, Amalthea or 
Thebe. This because the latter satellites have sizes of, at most, 100 km and masses which are 10~ 4 
- 10~ 6 the mass of Io. This choice takes into account the well known fact that the Titius-Bode rule 
works well for quite large and quite massive objects. For example, it does not work for comets or 
light asteroids. This criterium will be adopted also in the forthcoming considerations about the 
Saturn and Uranus systems. 



Saturn: 



2 e 2A Rin-ring = 2 ■ 1.5967 • 66 000 = 210 700 km 



Here we have used the radius of the ring "D", the inner one. We see that the number obtained 
lies half a way in between the orbits of Mimas (185 520 km) and Enceladus (238 020 km), which 
are the first two "big" satellites. The error is in the range of 11.5 % — 13.5 %. If we use the radius 
of the ring "C", namely Ri n - r ing = 74 500 km, we get for r\ = 237900 km, which agrees almost 
perfectly with the orbital radius of Enceladus (less than 1% error). However, Enceladus is not 
the first satellite (as regard the distance from Saturn) but only the second. Here also, as already 
done for Jupiter, we have discarded the too light bodies, and the satellites discovered only with 
spacecrafts (therefore, very small). The reasons are the same as in the above. 



Uranus: 



2 e 2A Rin-ring = 2 ■ 1.4662 • 41 840 = 122 690 km 



This value agrees with the orbital radius of Miranda (129 780 km), the first "big" moon consid- 
ered, within an error of 5.5%. Here we considered as first inner ring the ring "6". Miranda is the 
first satellite with relevant mass and size. In fact it was discovered from the Earth by Kuiper in 
1948. If we use the ring "Alpha" as inner ring, we get for the radius of the first satellite orbit 
r\ = 131 140 km, in agreement within 1 % with the Miranda orbital radius. 

Even if the agreement between the predicted inner radius of the rings and the observational 
data is not perfect (however with errors around 10 % or less), and the statistics of only three cases is 
really poor, nevertheless this "prediction" seems to corroborate the quantum- like model presented 
in this paper. On the other hand, it should be noted that these errors are of the same order of 
magnitude of those affecting the phenomenological Titius-Bode rule, therefore acceptable in this 
context. 
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11 Permitted orbits, dissipation and gravity 



One of the main objections that it is possible to rise against the existence of permitted discrete 
orbits in a planetary system is the following: It is a common experience, in this era of space travels, 
that a satellite can be put in any orbit we wish around the Earth, the Sun, or any other planet. Why 
therefore there should exist stable permitted orbits? In what sense they are "permitted"? How they 
are reached? Answers to these questions are suggested by the concepts of (chaotic) dissipation and 
limit cycles, emphasized in papers pQ (see also further references therein). Although in the proto- 
planetary nebula dust, particles and other bodies can be found at any distance from the central 
body, after a huge amount of time, friction, viscosity, and chaotic gravitational effects (for example 
planetary resonances) produce a dissipation of the initial state and bring matter to stabilize in 
several orbits, the limit cycles, where particles and dust aggregate to form planets. It is in this way 
that we can speak about "permitted orbits": they are the "limit cycles" of the dissipative processes 
started in the primitive nebula. 

Of course, it is in principle possible to put, today, a body in any orbit we wish. But if we wait for 
a time of the order of 5 • 10 9 years, and if such body has a sufficient mass, it is likely that we finally 
find it in one of the permitted limit cycles. 

The present model is therefore somehow complementary to the Classical Mechanics explanation that 
orbits in between the observed one are unstable, due to three-body forces (gravitational resonances). 
The quantum- like language of a wave equation is a tool to describe these collective effects of classical 
forces acting on unstable orbits. This "quantum" collective behaviour is an emerging phenomenon 
from an underlying classical deterministic (even if chaotic) dynamics. The description through a 
wave function is moreover very economical: it enables to get, at once, the satellite orbits and the 
rings. 

12 Conclusions 

In this paper we have shown that a wave equation is capable, under certain restrictive conditions, 
to describe some basic properties of the planetary systems, namely the law of the distances of 
the planets (or satellites) from the central body. The wave equation adopted for this scope is a 
deformation of the Schrodinger equation. 

We have been forced to choose a deformed Schrodinger equation by the analogy between the 
mechanism devised in [1] to produce quantization at atomic level via dissipation, and the chaotic 
dissipation occurred during the history of the proto-planetary nebulae. An dissipation mechanism 
could have been at work (of course, on much larger time scales) during the evolution of the planetary 
systems. From the primitive nebula, where all the orbits were filled by dust and rubble, we arrive, 
after 5 billion years of evolution, to the "quantized" orbits of today. 

Of course, having marked the analogies, also the evident differences must be underlined. We do 
not have quantum jumps in the Solar System, or quantum interferences between planets, neither 
quantum superpositions nor zero-point energy. A planetary system is not, of course, a quantum 
system. On the contrary, we have shown that a deformation of the Schrodinger equation (one of the 
basic tools of Quantum Mechanics) seems to be able to play a role also in the description of some 
quantum-like features of the planetary systems. Although indirectly, this quantum-like description 
seems also to corroborate the 't Hooft ideas on the origin of quantization from dissipation. 
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Appendix 

Linearity, self-adjoint ness and explicit form of = —isV v 
V v is defined on the orthonormal basis u m = e tm ^ /y/2n of £ 2 ([0, 2ir]) as 

P e im< ^ := i e mX e im ^ . 

Defining p^ = —isP v we have 

P(j) e im * = se mX e im * . 

Hence u m = e %m ^ j\f 7 hx are the eigenvectors of p^ with the eigenvalues \x m = se mX . 
Moreover we define Vtp to be linear by stipulating 

-p^^ e in<f> _|_ e im4>^ ._ ^ g nA e in<f> _|_ ^ e m\ e im4> _ -p^^inip^ _|_ •p^( e im 0) 

f^ae in ^) := ia e nX e in ^ = aV v {e in ^) . 

Being P^ linear on an orthonormal basis, then Vw is linear all over £ 2 ([0, 2tt]). 
We remind also the 

Theorem: If a linear operator is self-adjoint on an orthonormal basis of a Hilbert space H, then it 
is self-adjoint over all 7i. 

Therefore we have simply to show that p^ is self-adjoint on the basis u m . In fact we have 



2n se mX '"^ 7r 



h u m \u n ) = J # (p^u m y Un = — J #e- jm V^ = se mA 5,, 



and 



Jo Z7r Jo 



Besides, we check the identity 



In fact, since 



P v = ie 
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we have 



<""'• - { 1 • ( i\ Hi, + ... + l(-a) n d$ + ...^ • e im<i 

\n m n \ 

1 + Am + ... + + ...]• e im ^ = e mA e ^ 

n! 



Hence, the thesis. 

Confluent hypergeometric function 

The confluent hypergeometric function is defined via the series 

„, x a z a(a + 1) z 2 . 

^^^) = 1 + 7l! + WTl)2! + -' (12 - 8) 

which converges for all finite z; the parameter a is any number in C, and the parameter 7 must be 
different from zero and from any negative integer. If a is a negative integer (or zero), the function 
F(a,j,z) becomes a polynomial of degree \a\. 
The function F(at,'y,z) is a solution of the differential equation 

zu" + (7 - z)u -au = 0, (12.-7) 

as can be directly checked. 

With the substitution u = z x ~^u\^ eq. (|12.-7p transforms into 

zu" + (2 - 7 - z)u' x - (a - 7 + l)m = . 
From here we see that, for a non integer 7, eq. (|12.-7H admits also the integral 

z 1_7 F(a-7 + l, 2-7, z) 
which is linearly independent from (|12.-8p . so that the general solution of eq. ()12.-7p has the form 

u = c\F(a, 7, z) + C2 z 1 ~ J F(a — 7 + 1, 2 — 7, z) . 
The second term, contrary to the first, has a singularity in z = when 7 > 1. 
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